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The microscopic origin of the pseudogap state which exists in the underdoped cuprates remains 
unknown. The c-axis properties in the pseudogap regime are particularly anomalous. We use a 
recently proposed model of a d-density wave which leads to staggered currents and the doubling of 
the unit cell to investigate the c-axis kinetic energy and the optical sum rule. The density of states 
expected in the model is also considered. 



I. INTRODUCTION 



Pseudogap signatures in the underdoped regime of the. 
high-T c superconductors appear as intriguing featuresEl 
which could contain the key to a microscopic under- 
standing of the physics of this important region of the 
phase diagram. It is in this region that proximity to the 
Hubbard-Mott insulating phase would make one think 
that the underlying normal state is most profoundly dif- 
ferent from a usual metallic Fermi liquid state. Many as- 
sociated anomalous, properties are observed including im- 
portant violationcrB of the c-axis optical sum rule which 
related the missing area under the real part of the con- 
ductivity in the superconducting state, as compared to 
the normal state, to the corresponding c-axis superfluid 
density. Above the superconducting critical temperature 
T c , in the pseudogap regime, important changes in the 
plasma frequency are also observed. The preformed pair 
modejj which envisions that electron pairs are formed at 
the pseudogap transition temperature (T*) with T c the 
temperature at which phase coherence among the pairs is 
established, has helped understand some of these prop- 
erties. In particular the pseudogap phase corresponds 
to phase fluctuations with pseudogap energy the same 
as the pairing energy. An alternative to ihe preformed 
pair model elaborated upon by Chen et alB envisions ad- 
ditional incoherent pair excitations of finite momentum 
q which goes beyond the usual BCS formulation.-pf su- 
perconductivity. Very recently Chakravarty et alB have 
made a new proposal in which the pseudogap is associ- 
ated with the formation of a competing d-density wave 
(DDW) state on a nested Fermi surface (FS). This DDW 
order has important consequences such as staggered cur- 
rents with associated orbital magnetic moments which 
breaks parity and time-reversal symmetry. Since the 
original proposal, several works have appeared in which 
various properties associated with the DDW state have 
been elaborated upon with the aim of testing the vahcL 
ity of the model by comparing to experimental data.LTu 
ARPES data on the pseudogap have revealed it to have 
d-wave symmetry and to exist at least in the underdoped 
regime. The DDW order requires a nested FS and so 
should be most stable at half filling in the simplest of 



tight binding bands, and by assumption the DDW or- 
der is taken to have d-wave symmetry with maximum 
gap at (it, 0) and zero gap on the diagonals of the Cu0 2 
Brillouin zone . It is therefore of considerable interest 
to understand how much more of the known pseudogap 
physics such a model can explain. 

Charge transport along the c-axis has long been stud- 
ied in the cuprates and found to be anomalous. While for 
optimally doped YBa 2 Cu 3 07_5(YBCO) the c-axis DC 
resistivity tracks fairly well the linear in temperature (T) 



nderdoped case it 
with decreasing T 
Also the estimated 



in-plane resistivity, above T c , in th 
shows a semiconductor-like increai 
and has a maximum at T = T, 
mean free patbli3 becomes shorter than the interlayer dis- 
tance indication that Bloch transport is unlikely to ap- 
ply. A large decrease in the real part of the conductivity 
[<7i(u;)] is also observed at small frequency (uj) starting 
at the pseudogap temperature (T*). The energy scale for 
the pseudogap is of the same order as that for the super- 
conducting gap in YBa 2 Cu306.6- At T = 0, the suppres- 
sion of <j\ (id) for frequencies id below the gap scale is very 
well developed. In fact the additional spectral weight lost 
under the <7i(w) curve on entering the superconducting 
state is found to be smaller than the superfluid density 
determined at T ~ from the imaginary part of the con- 
ductivity. The c-axis sum rule is observed to Jae closer 
in value to a half than to its conventional valueEJ of one. 
This indicates that there must be an important change in 
kinetic energy on going from the pseudogap state at T c to 
the superconducting state at T = 0. The observed value 
of the c-axis sum rule of a half can be easily understood 
within a preformed pair model,li-3 as can the important 
changes in the optical plasma frequencies which are ob- 
served in the temperature range from T c to T* in the 
pseudogap state. Within the preformed pair model these 
effects are the result of the phase fluctuations of the su- 
perconducting order parameter. In this paper we study 
these issues within the assumption that the pseudogap 
state is due instead to the formation of the DDW order. 
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II. FORMALISM 

We begin, with a phenomenological model 
HamiltoniarEj of the form: H = —t^P. . CjtC 1(T — 

ME i(T n i* + J E<i,j> [ s i • s i - ; where C£ cre- 

ates a spin a electron on the site i, t the in- plane hopping 

amplitude, fi the chemical potential, and n, the occupa- 



tion on the site i. Si is the spin and J (> 0) a strengh 
of an exchange interaction. Based on the mean field 
approximation, defining A a = —J < diCi+s^ > and 
making an ansatz (— 1)*W<5 = —J < C^.Ci+i jCT >, we 
obtain the effective Hamiltonian for a combined super- 
conducting state with gap Ak and pseudogap state (with 
pseudogap Wk) 



Hmf = ^(e k - p)C+ a Ck* + [A k C+C- ki + h.c\ + ^ iW k C+ a C k+Qo 



(1) 



where £k = — 2t[cos(fc a; ) + cos(k y )], Ak = 2Ao[cos(fc a; ) — 
cos(fc y )], and Wk = 2Wo[cos(k x ) — cos(fc y )]. Here ek is 
the electron energy dispersion for a simple tight binding 
band with the FS at half filling coinciding with the anti- 
ferromagnetic Brillouin zone. The superconducting gap 
amplitude is Ao, and Wq is the corresponding amplitude 
for the pseudogap characterized by the commensurate 
wave vector Q = (tt,tt). For our simple band, therefore, 
£k+Q = -£k, A k+ Q = -A k , and Wk+Q = -Wk and 
this symmetry can be used to limit the summation on 
momentum to half the Brillouin zone. 

One can introduce an intersite interaction 
(V/2) E<i j> a a' n i<? n jo' to allow for a difference in in- 
teraction to exist in the d-wave superconducting (DSC) 
and DDW channel: J -> J - V for DSC and J -> J + V 
for DDW. Note that the site-dependent factor (— l) 1 in 
the ansatz made above ensures that nesting takes place 
for half filling and that Wg is pure-imaginary. 

In the 4x4 matrix form, Hmf = Ek ^k^kCk, 
where Ek means a sum over half of Brillouin zone, and 

C£ = ( C WV C - k l> C £+QV C -*-Ql)- Th e 4 x 4 matrix 
/ik is 

/ e k - n A k iW k 

A k -(eu-M) iW k 
iWk+Q e k +Q - pi A k +Q 
V iWk+Q A k+ Q -(e k +Q-ju) 

(2) 

This matrix defines our problem. We begin by studying 
c-axis properties and in particular the optical sum rule. 
To do this it is necessary to have some model for the 
c-axis charge transfer. For simplicity we will use 



(3) 



where t± is the transfer integral for electron hopping 
between two adjacent planes labeled by 1 and 2, re- 
spectively. We could equally well have taken incoher- 
ent coupling. The c-axis conductivity at frequency u) is 



<r c (0,w) = (t/uj)\u ret (0,u;)~e 2 d(H c ) 



an electron charge and an interlayer spacing, respectively, 
and 



11(0, iu) = 2e 2 t\dT Tr [^( k ' iu> n)G(k, iu) n 



and 



(H c ) = 2t\T J2 Tr [MG(k, iuj n )MG{k, 



(4) 



(5) 



In the above equations, the 4x4 matrix Green function 
- 1 and M - ( ^ ° 



f 3 



, where T3 



G(k,iu; n ) = [iu) n - hk] 

is a Pauli matrix. 

The c-axis superfluid density (or stiffness) (p s ) in the 
superconducting state is related to the limit as ui — > of 
the imaginary part of the conductivity. It is also related 
to the missing spectral weight under the real part of the 
conductivity on going from normal to superconducting 
state, which we denote by AW and define as 



AW • 



0+ 



duj [ff^M-ofM] 



(6) 



Here oj c is a cut-off frequency of order the bandwidth, 
and N and S stand for normal and superconducting state, 
respectively. Now, we have 



AW - 47re 2 d 



(H c ) s - (H c ) n 



= 4-7T lira [cjlma c (0,a;)] 



(7) 



From the above equations we obtain a sum rule in terms 
of Green functions G and Go = G(A = 0) as follows: 



aw 1 , !E:, k E,., g§,«] 



Ew.k Ei,j PijGij 



where e and d are where 
I 



^ a ij Go,ij] — (^11 — ^0,ll) — (^13 _ ^0,13) - (^24 — ^0,24) + (^33 - ^0,33) 



(8) 



(9) 
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and 



E* 



r 2 



G? 9 



a 2 



G 



23 



G 



34 



Note that the superfluid density p s is given by 



(10) 



(11) 



where C = 327re 2 i 2 L <i. The equations above are easily 
generalized for incoherent interlayer coupling. As T — » 
(zero temperature limit), 



J 



C E A k 



£ 3 



£ 3 

^2k 



(12) 



w here g lk = ^(^ok - W + A k 
V(^ok + /i) 2 + A 2 with S k 



and E 



2k — 



|£k|- 

-^2k 



W 2 . 

In the case when the DDW order W^ = 
For 



0, R 



Ok 



e k > 0, E n 



and 



(e k+Q - pf 



A k+Q while for e k < Ei k 
and E 2 k 



(ek+Q - ^) 2 + A k+Q and E^k = \/(ek - m) 5 
Here we have used the properties associated with the 
commcmsurate wave vector Q = (tt, n). Then, we obtain 



E 

k 



A k 



A 2 

^k+Q 



[(e k -/i) 2 + A 2 ] 



2l3/2 



(e k +Q - l±f 



A 2 

^k+Q 



3/2 



^E 



A| 

K 



(13) 



where the summation over k is unrestricted i.e. is over 
the entire Brillouin zone. This is the usual expression of 
p s with Ek — y/ (e k — p) 2 + A k for the pure DSC case 
with no DDW order. 



III. C-AXIS RESPONSE 

Based on experimental observations!^ we choose the 
doping dependence of the DDW gap amplitude Wq(x) 
to be Wq(x) = 0.04(1 — x/x c ) in unit of t. Here x 
is doping and its critical value is x c — 0.2, where the 
DDW gap disappears. Also we assume Ao(x) = 0.02 for 
0.05 < x < 0.2. A small variation of A (x) near the opti-. 
mal doping (x ~ 0.16), which is seen in the experiment ,t£l 
is not important because our main consideration is fo- 
cused on the underdoped regime of the cuprates. [See a 
later discussion and an inset of Fig. 2.] In Fig. 1, we plot 
the calculated c-axis superfluid density vs the amplitude 
of the pseudogap Wq(x) at T = assuming a constant 
chemical potential. Such an assumption is of course not 
correct since the chemical potential also changes with 
the dopping x; however, keeping it fixed will help us to 
understand the physics. The solid curve is the c-axis su- 
perfluid density with the chemical potential p = —0.06. 
The dashed curve is for a larger value of /i = —0.2. OD 
and UD at the top of Fig. 1 stand for overdoped and un- 
derdoped regime, respectively. One can see from Fig. 1 
that when the absolute value of the chemical potential 
is smaller than the pseudogap value of W(x) — AWq(x), 
which is its maximum for a given doping x, then the su- 
perfluid density is reduced as W{x) in increased. This 
happens to be the case for most of the doping values x 
between 0.05 and 0.2 in the solid curve. On the other 
hand, for larger values of \p\ the opposite holds (dashed 
curve). To understand the physics behind this behavior 
we return to the well known result embodied in Eq. fllE 



r 



which gives the superfluid density p s o for a pure d-wave 
superconductor with no DDW gap i.e. Wo(x) — 0. In 
this case assuming a cylindrical FS p s o = CN(0)/2, where 
N(0) is the density of states at the FS. In this simple 
model, N(0) is constant throughout the band; however, 
if it did vary, then it would be some appropriate energy 
average of 2V(u>) around the FS on the scale of A that 
would replace it. We stress that for p s0 the size of the zero 
temperature superfluid density is independent of the size 
of the superconducting gap Ao and depends instead on 
normal state parameters. While Eq. (|2|) has an explicit 
factor of A k in the numerator, this does not translate 
directly into a reduced superfluid density at T = as 
Ao(a;) decreases. 

As the pseudogap develops it competes with the super- 
conductivity but we can understand the trends seen in 
Fig. 1 for p s as a function of Wq(x) by considering, first, 
the effect of Wq (x) on the electronic density of states and 
then thinking of switching on the superconductivity. In 
the inset of Fig. 1, we show the quasiparticle density of 
states for the pure pseudogap state (i.e. A (x) = 0). 
The solid curve is for Wo = 0.03 i.e. x = 0.05 and the 
dotted curve is for Wo = 0.016 i.e. x = 0.12. We will re- 
turn to a detail discussion of the density of states (DOS) 
in a later section. For now suffice it to note that the 
logarithmic si ngularitie s associated with the DOS are at 
\n\ ± AWot/yJt 2 + Wo 1 . For a wide band with t » W 
they would be at \p\ ± 4Wq. Now for fi = —0.2 the 
Fermi surface FSi is localted away from the DDW gap 
and the DOS aound the Fermi energy increases as dop- 
ing ranges from overdoped to underdoped i.e. as Wq(x) 
increases. Consequently, p s at T = increase. However, 
when p = —0.06 the Fermi surface FS2 sits inside the 
DDW gap shown for most values of x and an average of 
the DOS around FS2 is reduced as x becomes smaller i.e. 
Wo (a;) becomes larger. Now the value of the superfluid 
density decreases as the underdoped regime is entered. 
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Thus the opening of the pseudogap will decrease p s at 
T = only if \p\ is small enough in the DDW model. For 
realistic values of x, however, the chemical potential p is 
not close to zero when it is calculated self-consistently for 
the simple tight binding band structure we have assumed 
and in the presence of a DDW order. While one can take 
empirical expression for the doping dependence of super- 
conducting gap and pseudogap , an arbitray choice of 
fx(x) would not be consistent with our model band struc- 
ture and, therefore, is not allowed. We must calculate p 
from the filling. 

From this discussion we conclude that the opening-up 
of a pseudogap is not sufficient to lead to a large re- 
duction in the c-axis superfluid density. On the other 
hand the interlayer hopping matrix element t± (x) in the 
YBCO series is known to decrease almost exponentially 
as the doping is decreased. This factor will dominate over 
effects of the DDW order in consideration of the suppres- 
sion of the c-axis superfluid density with doping. These 
considerations imply that the doping dependence of the 
c-axis superfluid density is not a good quantity in which 
to study the role of the DDW order. A better choice 
is the c-axis optical conductivity sum rule because it is 
independent of the magnitude of the interlayer hopping 
amplitude and, therefore, of its dependence on doping. 
We next turn to the calculation of the c-axis sum rule and 
focus on the issue of whether or not the DDW model can 
describe the experimental observations .□ In the YBCO 
system a conventional sum rulalj of one is observed for 
the optimally doped case and of about a half for an un- 
der doped sample. 

For the sum rule calculation we will make use of some 
experimental observation rather than proceeding to a 
complete self-consistent calculation. This is reasonable 
since the nature of the interaction J — V and J + V, 
which determine the size of the amplitude of DSC and 
DDW gap respectively in the formalism, is not known. 
Assumptions on the variation of Ao and Wq as func- 
tions of x based on empirical expressions correspond to 
specific assumption about the unknown variation of the 
above interaction parameters. However, in order to take 
into account the band structure and doping we deter- 
mine the chemical potential as a function of doping and 
temperature by solving the filling equation derived in the 
Bogoliubov-de Gennes formalism. 

In practice a more useful expression than Eq. (|8|) to 
calculate the sum rule, or the normalized missing spec- 
tral weight (NMSW) AN/p s is 



J 



n = 1 



It, 



£ 2k V 2T 



AN 4ire 2 d 
= 1 + 

Ps Ps 



[H c y - (H c ) n 



(14) 



Of course Eq. (g) is more direct if one wishes to under- 
stand how a sum rule of a half can be obtained. This 
results when the Green's functions in the second term 
of the right hand side cancel between superconducting 
and pseudogap state as in the preform pair model.E-3 We 
see from Eq. ([IJ) that the kinetic energy difference be- 
tween the superconducting and pseudogap state divided 
by the superfluid density determines the c-axis conduc- 
tivity sum rule. Mathematical expressions for the kinetic 
energy and the superfluid density are as follows: 



and 



4ne 2 d{H c ) = 



CT EE 



7i +72 
P 2 



if 2 



(15) 



(16) 



where 



7i = K - E < 



Ok 



„ 2 )(to 2 n + E 2 k + Al + p 2 y (17) 



72 = V^ k (£& + 3o£ + 3A£ + p, ) 
(3 = (w* + E 2 k + A£ + p 2 ) 2 — 4p 2 E 2 k 



Ok 



Al + p 2 Y + 4p 2 E l 



2^2 
Ok 



(18) 
(19) 
(20) 



The kinetic energy of each state should presumably be 
calculated at zero temperature because in some notion 
of the theory of interlayer coupling, the kinetic energy 
difference is associated with the condensation energy of 
superconductors.tll However, Basov et al. have taken 
T ~ T c since it is not easily to access the pseudogap 
state at low T. In the calculation of the kinetic energy 
we choose T = 0.01 for the superconducting state and 
T = 0.1 for the pseudogap state. If t = 2000JT, then the 
working temperatures are 20K and 20QK in the super- 
conducting and pseudogap state, respectively. 

As we mentioned ealier, we determine self-consistently 
the chemical potential p for a given doping from the band 
structure as well as the given temperature even though 
we adopt the experimental observation for Wq and Ao. 
The expression for the filling n, in the self-consistent for- 
malism, is written as: 



^L^tanhf^l 



(21) 



In Fig. 2, we plot NMSW (AN/p s ) as a function of dop- 
ing (x) . In the inset of Fig. 2 we reproduce some of the 
experimental phase diagram obtained recently by Tal- 
lon and Loramtj which shows critical temperature (plus 



r 



sign +), pseudogap value Wo (solid square) and super- 
condicting gap (solid triangle) as functions of doping for 
Yo.8Cao.2Ba2Cu307_,5. Here the details are not impor- 
tant and we take the data as typical of the cuprates. We 



4 



approximate, in unit of t, Wo(x) — 0.04(1 — x/x c ) and 
A (x) = 0.02 at T = 0.01, and W {x) = 0.025(1 - x/x c ) 
with Aq(x) = at T = 0.1. As shown, in the underdoped 
regime NMSW is less than one while it saturates to about 
one beyond the optimal doping. Different values of Wb(0) 
at T = 0.1 do not change the overall behavior of NMSW 
vs x. However, we did find that a different choice of 
working temperature for the pseudogap state can change 
the doping dependence of NMSW and, thus, the behav- 
ior shown in Fig. 2 is not robust. It has been found that 
a temperature dependence of NMSW for a d-wave su- 
perconductor with a cylindrical FS was negligible. The 
behaviors of NMSW in the tight binding band with a 
pseudogap are different. While the DDW model gives 
the right trend for the doping dependence of the c-axis 
sum rule it does not fall to the value oLa half observed in 
the experimental work of Basov et al.a We did not find 
a set of parameter with which this could be achieved. 
In this regard we should mention that in the preformed 
pair model it is argued that the result AJV/p s = 1/2 
in the underdoped regimejij follows simply because the 
Green's functions appearing in the numerator of the sec- 



J 



ond term in Eq. (g) (which are diagonal in this case) are 
the same for superconducting and pseudogap state. In 
the preformed pair model only the off-diagonal part of 
the Green's functions differ in the two regime. It is this 
part which deals with fluctuations of the phase coherence 
of the Cooper pairs. 

Next we turn our attention to the pseudogap state 
above T c but below the pseudogap temperature T* . In 
this regime the optical spectral weight or plasma, fre- 
quency is observed to vary with temperatureBEjEj We 
reproduce in the inset of Fig. 4 the data on YBa2Cu306.6 
given in Ref.llil This implies that a change in c-axis ki- 
netic energy occurs as temperature is lowered from T*, 
and in our model as the DDW gap grows. 

Within a self-consistent formalism, we calculate the ki- 
netic energy above T c and below T* in the regime when 
DDW gap is included in the calculation, and we compare 
the theoretical calculations with experimental results.E-3 
For simplicity, we do not consider inhomogeneity in the 
system. The self-consistent equations for Wu and the 
filling n are: 



k 



Worj 



tanh 



f #0k + /i 

I 2T 



tanh 



f ^ok - M 
I 2T 



(22) 



and Eq. (21) with A k = 0, where Vddw is the DDW 



channel interaction and 7/ k = cos(fc x ) — cos(k y ). For 
given values of Vddw and n, Wo(T) and fJ.(T) are self- 
consistently determined. In the calculation we choose the 
DDW channel interaction Vddw — 0.945 and the values 
n = 1, 0.99, 0.97, and 0.95 for the filling. In fig. 3, we plot 
Wo(T) versus T in unit of t. In the inset we show n(T) as 
a function of T for n = 0.95. In this case, if we assume 
t = 2000 AT, then T* ~ 215 AT. Since below T r , which 



J 



{Ho) 



k 



is taken as 80A, we would have to include an equation 
for Ak as well in the self-consistent calculation and we 
have not, we simply indicate Wq(T) as a dashed line for 
T < T c . Using W (T) and fi(T) for different n's, we cal- 
culate the normalized kinetic energy (H C )(T)/(H C )(T*) 
for T c < T < T*. An analytic expression for (H c ) is 
achievable only if the magnitude of the chemical potel- 
tial \n\ is much less than a working temperature; namely, 
\n\/T <C 1. In this instance it can be shown that 



tanh 



V T 



3 tanh 2 (^)-l 



(23) 



r 



The correction to the above equation is of order (/i/T) 4 . 
In Fig. 4, we plot our numerical results for the normalized 
kinetic energy for T c < T < T* . The curves are labeled 
by the value of filling n with n = 1.0 half filling and 
n = 0.95 corresponding to a value of doping x = 0.05. 
It is clear from the figure that near half filling the c-axis 
kinetic energy and so the optical spectral weight does de- 
crease very significantly as we go from T* toward T c . The 
reduction for half filling is by a factor of 5 more than nec- 
essary to agree with experiment for YBa 2 Cu 3 06.6 shown 
in the inset of Fig. 4 which gives the experimental results, 
for the square of plasma frequency as a function of T.til 
However, as we move away from half filling, the initial re- 
duction with decreasing temperature out of T* is rapidly 
suppressed. Also at the lower temperatures considered, 



the optical spectral weight begins to increase again. The 
tight binding band structure we have used in our calcula- 
tions is at best genetic for the oxides and may not apply 
in a quantitative sense for a particular case. Nevertheless 
it is clear that the model of DDW order as an explana- 
tion of the physics of the pseudogap regime gives results 
for the change in c-axis kinetic energy with temperature 
which are rather sensitive to the assumed parameters, for 
example to doping. 



IV. DENSITY OF STATES 

Renner et aL@ have studied the evolution of the pseu- 
dogap features in STM tunneling (SIN) in a series of 
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I^S^CaCu^Os+a (Bi2212) as a function of doping from 
underdoped to overdoped regime of the phase diagram. 
Somewhat complimentary data are given in DeWilde et 
al.tl and in Miyakawa et allEn where the authors con- 
centrate more on the hump and dip feature seen around 
twice the gap energy in SIN and, at three times in Joseph- 
son junctions. Renner et alHl find a superconducting 
gap which is nearly temperature independent up to T c 
at which point it merges smoothly into a second gap like 
feature centered at the FS. They find the pseudogap to 
be present both in underdoped and overdoped samples 
and that its size scales in magnitude with the super- 
conducting gap. This argues for a common origin and 
these authors conclude that the data is consistent with 
the idea of preformed pairs. Because they find that the 
pseudogap is tide to the FS the data do not support a 
conventional band structure explanation. On the Qthps 
hand recent intrinsic tunneling spectroscopic resultsE3'Efl 
for mesa on Bi2212 single crystals (SIS) have revealed 
distinct features that can be associated with the super- 
conducting gap and with the pseudogap. The tempera- 
ture dependence of the superconducting peak structure 
in the dynamic conductance shows that it closes at T c 
while the hump structure which is associated with the 
pseudogap is unaffected by the superconducting transi- 
tion. This was taken as evidence that the two phenomena 
are distinct and that the data do not favor a preformed 



pair model. Later magnetic field studies showed that the 
structure identified with the pseudogap is insensitive to 
magnetic field while the superconducting gap is strongly 
suppressed by the field,Ejc3 which is taken as further 
confirmation of the identification made. A feature of the 
data which is relevant to our work is that the supercon- 
ducting gap corresponds to sharp peaks which are seen 
to grow even sharper with decreasing temperature and 
fall inside the pseudogap humps. The DDW model does 
not relate gap and pseudogap directly and so is consis- 
tent with the above interpretation. In the self-consistent 
phenomenological model the two are related to J — V 
and J + V respectively and, therefore, they can be quite 
different in size and in variation with temperature. On 
the other hand, the model does not give superconduct- 
ing gap falling inside a larger pseudogap structure. This 
behavior is generic to the model and has its origin in the 
fact that the superconducting gap opens up at the FS 
while the DDW gap is centered at the antiferromagnetic 
Brillouin zone boundary as described below. 

Although the density of states (DOS) N(uj) is related 
to the in-plane dynamics, it is nevertheless relevant to 
our present discussion because we have already seen that 
the c-axis response reflect in-plane dynamics. Employing 
standard manipulations, it can be shown that the quasi- 
particle DOS N(u) is: 



iYH = -2^Im[Gf 1 + G|] 

k 

= 2 ~ E lk ) + v 2 lk 5(uj + E lk ) + u 2 J(uj - E 2k ) + v^S(uj + E 2k )] 



(24) 



Here a factor of 2 is for the summation over spin, 
G R is the retarded Green's function and the coher- 
ent factors are uf k = | [1 + (E ok — A*)/ ^ik] and u 2 k = 
\ [1 - (£ok + M)/£ 2 k] with i4+t& = 1 (i = 1, 2). There 
are three limits of special interest: 

i) When W k = 0, one can show that N(u>) reduces 
to DOS for the pure DSC case as we already showed 



for the superfluid density. In this instance, for e k > 0, 



J 



-Elk 
u lk : 



k, ^2k = E k+Q , u{ k = i[l + £k/#k], and 
| [1 + Ck+Q/^k+Q], where £ k = e k - M and E k = 
If e k < 0, then # lk = E k+Q , E 2k = E kl 



'ik 



= I ri 



■£k+Q/£k+Q], and u 2 k = i[l + £k/£] 



Therefore, we obtain 



N{cu) = £{[1 + W^k] S(u - E k ) + [1 - £k/£k] 5{lu + E k ) 

k 

+ [1 + Ck+Q/Sk+q] S(ui - E k+Q ) + [1 - &+Q/£k+q] S(u + £ k+Q )} 
= £{[1 + W£k] 8{u - E k ) + [I - £ k /£ k ] 6(lu + E k )}. 



(25) 



ii) When A k = 0, E lk = |JSok— A*| and E 2k = \E ok +fx\. N(ui) = ^ [S (w - {E ok - (j,)) + 6 (w + {E ok + (J,))] 
le coherent factors are now u\ k = | [1 + sgn(i? k — A*)] k 



The 

and u 2i( . = i [1 — sgn(i?ok + A 4 )]- Simple algebra shows 
that the DOS for a pure DDW case is 



(26) 



G 



iii) When p = 0, but, with nonzero DSC and DDW 
gap, Eik = E 2 k = S k , where S k = ^/e£ + A£ with 
Ak = \/A k + W£. In this case the DOS becomes 

N(u) = l S ( u ~ S k) + S(lu + Sk)] • (27) 

k 

Note that this DOS is nothing but the DOS of a pure 
DSC case or the DOS of a pure DDW case with a gap 
A k for (j, = 0. 

For the solid curve in the top frame of Fig. 5, Wq = 0.04 
and Ao = 0.02 with the chemical potential p = —0.25 
which is sufficiently large in absolute value that super- 
conducting gap and pseudogap are well separated, and 
there is distinct and characteristic structure associated 
with each. The superconducting gap froms at the FS 
(u> = 0) in the figure and the DDW gap at ui = \p\. The 
position in energy u) of the singularities corresponding 
to 4A and 4Wo are renormalized because both gap and 
pseudogap are present and interfere, and because p is not 
zero. But these renormalizations are not large. For the 
pseudogap alone (A = case) there should be peaks of 
equal heights at cu ~ \p\ ± 4Wq i.e. oj ~ 0.41 and 0.09. 
While the uppermost peak is shifted upward only very 
slightly by the presence of A to 0.418, the second peak 
is shifted more importantly to 0.12. If the pseudogap 
were zero, the superconducting peaks would be close but 
not quite at u> = ±0.08 because p is not zero in this exam- 
ple. They are at w = ±0.076. Other than this small shift 
the gap and pseudogap structures are quite separate in 
the solid curve. Note also from the general mathematical 
form Eq. (|24|) for the quasiparticle DOS that structures 
will fall symmetrically in energy about w = but they 
are of different height. In particular we note the small 
but clearly visible peak at oj = —0.12. The dot-dashed 
curve in the top frame of Fig. 5 is for comparison with 
the solid curve and differs from it only through a different 
value of Wq = 0.03 instead of 0.04. This change clearly 
shifts the the two prominant pseudogap peaks leaving 
the superconducting gap structure much less affected. It 
is clear from this comparison that gap and pseudogap 
are pretty independent of each other although each is 
somewhat affected by the presence of the other. If the 
chemical potential is reduced to p = —0.08 as in the 
middle frame of Fig. 5, we can see that the competition 
between gap anf pseudogap becomes much more severe 
particularly in the region around the FS and the chemical 
potential which are now closer together. The position of 
the most prominant upper peak falls at to = 0.173 above 
the chemical potential (\p\ = 0.08). This value is still 
close to 4Wq = 0.16. The third peak, which in the top 
frame was identified with the gap peak, falls at u> = 0.036 
which is a factor of 2 lower than 4Ao = 0.08 so that if 
we should identify this with the gap, we would have to 
conclude that it is strongly suppressed by the presence 
of the DDW order. This is expected since both DSC and 
DDW order compete for the available states in this re- 
gion of energy. This interpretaion finds further support 



in the fact that a change in the value of gap Ao affects 
the position of this third peak strongly while a change of 
the DDW gap Wq does not. The position of the second 
peak however changes significantly with a change in Ao 
or Wo so that this peak is a true mixture of both gap and 
pseudogap . 

In the bottom frame of Fig. 5, we show results for 
/i = i.e. at half filling. In this case, as we have seen 
in the D OS, gap a nd pseudogap become a single entity 
A = yAo + Wo = 0.045. A simple gap structure is 
obtained symmetric about oj = 0. In this case, the singu- 
larities are located at ±4Ao£/ *Jt 2 + Aq. It is clear from 
the trends exhibited in Fig. 5 that the DDW model does 
not allow for a separate well identified set of two gaps 
with superconducting gap falling inside a larger pseudo- 
gap about u> = 0. This could only happen if gap and 
pseudogap where both pinned to the FS. In the DDW 
model only the superconducting gap is pinned to the FS 
while the DDW gap opens up at the energy of the anti- 
ferromagnetic Brillouin zone. Of course the intrinsic tun- 
neling experiments do not measure the DOS directly as 
they involve SIS junctions. Our calculations can be coi&i 
pared more directly to the SIN data of Renner et a/.Ej 
which do not presently show the evolution with doping 
expected in a DDW model. 



V. CONCLUSIONS 

As a possible candidate model for the pseudogap state 
seen in the underdoped regime of the high-T c cuprates, 
we have studied the effects of the formation of Z?-density 
wave (DDW) order on several properties. While it is 
found that the value of zero temperature superfluid den- 
sity (p s ) is impacted by the opening- up of a DDW gap 
at the antiferromagnetic Brillouin zone, its effect can in- 
crease as well as decrease the value of p s . If the chemical 
potential is small and near zero, the effective normal state 
density of electronic states around the Fermi surface (FS) 
is on average reduced when the pseudogap is increase by 
going to a more underdoped case. This leads to a re- 
duction in p s . But if the chemical potential falls well 
above the pseudogap energy, the opposite holds. In any 
case these effects are not dominant and are masked by 
the near exponential reduction in interlayer transfer ma- 
trix element known to exist in YBa 2 Cu 3 07_,5 as one goes 
from overdoped (where t± ~ 30 — 40mev) tojinderdoped 
case (where t± can be a fraction of a meV).E£!l 

A quantity which does not depend explicitly on t± and 
its variation with doping is the c-axis optical sum rule, 
which is known to be anomalous in underdoped cuprates 
while it takes on its conventional value of one in optimally 
doped YBa2Cu307_,5. In this case we find that opening 
a DDW gap will affect the value of the ratio of the miss- 
ing area under the difference between the real part of the 
optical conductivity in pseudogap and superconducting 
state. When it is normalized to the value of the c-axis 
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superfiuid density, which is determined from the zero fre- 
quency limit of the imaginary part of the conductivity, 
it is found to be close to one in the overdoped regime 
and decrease towards a half in the underdoped case. In 
the case presented, the minimum value obtained was a 
little less than 0.8 but the calculations show that the size 
of the reduction depends on details of band structure, on 
filling, and on size of the gap and pseudogap as well as on 
the reference temperature used. A reduction can occur 
but its size is not a robust feature of the model. 

We also considered how the area under the real part 
of the conductivity which gives the plasma frequency (or 
optical spectral weight) is affected by the introduction of 
the DDW order. It is found that near half filling, the op- 
tical spectral weight is substantially reduced as the tem- 
perature (T) is decreased from the pseudogap tempera- 
ture (T*) towards the superconducting transition tem- 
perature (T c ). Agreement with experiment is possible 
near half filling. When the doping is increased, however, 
away from half filling, the plasma frequency is found to 
have a minimum in the region T c < T < T* . This be- 
havior is different from that expected on the bases of 
the preformed pair model where the phase fluctuations 
are the cause of the variation in optical spectral weight 
which only decreases with decreasing T. 

To acheive some physical understanding of our c-axis 
result we found it useful to introduce and consider the 
modifications brought to the in-plane quasiparticle den- 
sity of states through the growth of DDW order. Be- 
cause a generic feature of the model is that superconduct- 
ing gap opens at the Fermi surface while the DDW gap 
opens instead at the magnetic Brillouin zone, the result- 
ing density of states N(oj) never shows superconducting 
gap features distinct from pseudogap feature which are 
also positioned inside the pseudogap humps as reported 
in a recent experiment. The prediction is that when the 
chemical potential (fx) is large enough, gap and pseudo- 
gap features are separated by \fx\ and are quite distinct. 
When fi = 0, gap and pseudogap can no longer be dis- 
tinguished and the square root of the sum of the squares 
plays the role of a single gap instead, reminiscent of the 
preformed pair model. 
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FIG. 1. The c-axis superfluid density p s divided by its maximum value as a function of the value of the pseudogap amplitude 
Wq(x) which depends on the doping x. In the solid curve the chemical potential \i is fixed at —0.06 and in the dashed curve 
at —0.2 in a unit of t. The inset gives the density of quasiparticle states (DOS) in a pure DDW state. The solid curve is for 
Wo = 0.03 and the dotted for 0.016. Two possible position for the Fermi surface (FSi and FS2) are indicated by vertical lines. 
OD and UD on the top stand for overdoped regime and underdoped regime, respectively. 

FIG. 2. The value of the c-axis sum rule AJ\f/p s as a function of doping x. The changes reflect the opening oJ-jthe DDW 
gap as the underdoped regimed is entered. The inset shows a typical phase diagram for the cuprates from Ref.Ej with the 
superconducting dome (plus sign +), pseudogap (solid square) and superconducting gap (solid triangle) in unit of [K]. Note 
that for both superconducting gap and pseudogap their maximum on the FS is four times Ao and Wo, respectively. 

FIG. 3. The value of the DDW gap Wo(T) as a function of temperature (T) determined self-consistently for a given value 
of parameters (n = 0.95 and Vddw — 0.945) as described in the text. The inset shows the corresponding chemical potential fj, 
as a function of T. 

FIG. 4. The c-axis kinetic divided by its value at T* as a function of temperature (T) in the range T c < T < T* . The curves 
are labeled by the values of the filling n and related to the doping p-|= 1 — n. In the inset we show the square of the normalized 
plasma frequency as a function of temperature T obtained in Ref.EJ for YBa2Cu306.6 

FIG. 5. The density of states (DOS) N(ui) as a function of energy u including both DDW order and superconductivity. For 
the solid curve in each frame, the DDW gap Wo = 0.04 and the superconducting gap Ao = 0.02 while the dot-dashed curve in 
the top frame has Wo = 0.03 instead. The frames differ in choice of chemical potential fj,. In the top frame fj, — —0.25, middle 
fj, = —0.08, and in the bottom fi — 0. The dahed line indicates to — \y\ and FS stands for the Fermi surface. 
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Fig.2 (Kim et al) 
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